home *** CD-ROM | disk | FTP | other *** search
/ QRZ! Ham Radio 8 / QRZ Ham Radio Callsign Database - Volume 8.iso / pc / files / dsp / 96000tar.z / 96000tar / 96000 / appb / b143a.asm < prev    next >
Assembly Source File  |  1992-04-28  |  5KB  |  97 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.43    Line Drawing  
  8. ; B.1.43.1    Floating-Point Incremental Line Drawing Algorithm  
  9. ;This algorithm generates points along a line given the endpoints. As the  coordinate along one axis is 
  10. ;incremented in fixed point, the other  coordinate is incremented in floating-point and then converted 
  11. ;to fixed  point.  A full line drawing algorithm which draws lines in all  directions is given below.  
  12. ;Registers: 
  13. ; 
  14. ;  d0 = temporary      d4 = temporary x1     d8 = 
  15. ;  d1 = temporary      d5 = temporary y1     d9 = 2.0 
  16. ;  d2 = x1 (dx)        d6 = x0 and xScreen 
  17. ;  d3 = y1 (dy)        d7 = y0 and yScreen 
  18. ; 
  19. ;                                                           Program    ICycles 
  20.                                                              Words 
  21. ; Calculate dx and dy 
  22.   fsub.s d6,d2  d2.s,d4.s                                ;  1    1 
  23.   fsub.s d7,d3  d3.s,d5.s                                ;  1    1 
  24. ; Determine whether to increment x or y 
  25.   fcmpm  d3,d2                                           ;  1    1 
  26.   fjge   _inc_x                                          ;  2    2 
  27. ; Switch endpoints if necessary 
  28. _inc_y 
  29.   ftst   d3                    d2.s,d0.s                 ;  1    1 
  30.   ftfr.s d4,d6                 fflt                      ;  1    1 
  31.   ftfr.s d5,d7                 fflt                      ;  1    1 
  32. ; Fix y0 and dy 
  33.   int    d7                    d3.s,d1.s                ;   1    1 
  34.   int    d1                                              ;  1    1 
  35.   neg    d1                    iflt                      ;  1    1 
  36.   jeq    _draw1_y                                        ;  2    2 
  37. ; Calculate dx/dy 
  38.   fseedd d3,d4                                           ;  1    1 
  39.   fmpy.s d3,d4,d5              d9.s,d2.s                 ;  1    1 
  40.   fmpy   d0,d4,d0 fsub.s d5,d2 d2.s,d3.s                 ;  1    1 
  41.   fmpy.s d5,d2,d5              d2.s,d4.s                 ;  1    1 
  42.   fmpy   d0,d4,d0 fsub.s d5,d3                           ;  1    1 
  43.   fmpy.s d0,d3,d0              d6.s,d2.s                 ;  1    1 
  44. ; Draw first point                                       ;
  45.   int    d6                                              ;  1    1 
  46.   jsr    _draw_point                        ;  application dependent 
  47. ; d0 = dx/dy  d1 = dy  d6 = x0  d7 = y0 
  48.   do     d1.l,_end_y                                     ;  2    3 
  49.   fadd.x d0,d2                                           ;  1    1 
  50.   inc    d7                    d2.s,d6.s                 ;  1    1 
  51.   int    d6                                              ;  1    1 
  52.   jsr    _draw_point                        ;  application dependent 
  53. _end_y 
  54.   rts                                                    ;  2    2 
  55.  
  56. _draw1_y 
  57.   int d6                                                 ;  1    1 
  58.   jsr    _draw_point                        ;  application dependent 
  59.   rts                                                    ;  2    2 
  60. ; Switch endpoints if necessary 
  61. _inc_x 
  62.   ftst   d2                    d3.s,d0.s                 ;  1    1 
  63.   ftfr.s d4,d6                 fflt                      ;  1    1 
  64.   ftfr.s d5,d7                 fflt                      ;  1    1 
  65. ; Fix x0 and dx                                          ;
  66.   int    d6                    d2.s,d1.s                 ;  1    1 
  67.   int    d1                                              ;  1    1 
  68.   neg    d1                    iflt                      ;  1    1 
  69.   jeq                          _draw1_x                  ;  2    2 
  70. ; Calculate dy/dx                                        ;
  71.   fseedd d2,d4                                           ;  1    1 
  72.   fmpy.s d2,d4,d5              d9.s,d2.s                 ;  1    1 
  73.   fmpy   d0,d4,d0 fsub.s d5,d2 d2.s,d3.s                 ;  1    1 
  74.   fmpy.s d5,d2,d5              d2.s,d4.s                 ;  1    1 
  75.   fmpy   d0,d4,d0 fsub.s d5,d3                           ;  1    1 
  76.   fmpy.s d0,d3,d0              d7.s,d2.s                 ;  1    1 
  77. ; Draw first point                                       ;
  78.   int    d7                                              ;  1    1 
  79.   jsr    _draw_point                        ;  application dependent 
  80. ; d0 = dy/dx  d1 = dx  d6 = x0  d7 = y0 
  81.   do     d1.l,_end_x                                     ;  2    3 
  82.   fadd.x d0,d2                                           ;  1    1 
  83.   inc    d6                    d2.s,d7.s                 ;  1    1 
  84.   int    d7                                              ;  1    1 
  85.   jsr    _draw_point                        ;  application dependent 
  86. _end_x                                                   ;
  87.   rts                                                    ;  2    2 
  88. _draw1_x                                                 ;
  89.   int d7                                                 ;  1    1 
  90.   jsr    _draw_point                        ;  application dependent 
  91.   rts                                                    ;  2    2 
  92.  
  93.  
  94. ;Performance: 
  95. ;  Trivial case: (single point)  16      cycles 
  96. ;  Other cases:                  25 + 3n cycles 
  97.